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The K-ban6 intensity profile of R Leonis probed by VLTI/VINCI 
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Abstract. We present near-infrared ^-band interferometric measurements of the Mira star R Leonis obtained in April 2001 and 
January 2002 with the VUTl, the commissioning instrument VINCI, and the two test siderostats. These epochs correspond to 
near-maximum stellar variability phases ~0.08 and ~ 1.02 (one cycle later), respectively. The April 2001 data cover a range 
of spatial frequencies (31-35 cycles/arcsecond) within the first lobe of the visibility function. These measurements indicate 
a center-to-limb intensity variation (CLV) that is clearly different from a uniform disk (UD) intensity profile. We show that 
these measured visibility values are consistent with predictions from recent self-excited dynamic Mira model atmospheres that 
include molecular shells close to continuum-forming layers. We derive high-precision Rosseland diameters of 28.5 ± 0.4 mas 
and 26.2 ± 0.8 mas for the April 2001 and January 2002 data, respectively. Together with literature estimates of the distance 
and the bolometric flux, these values correspond to linear radii of 350^^5} R© and 320^^5} R©, and to effective temperatures of 
2930 ± 270 K and 3080 ± 3 10 K, respectively. 

Key words. Techniques: interferometric - Stars: late-type - Stars: AGB and post-AGB - Stars: atmospheres - Stars: funda- 
mental parameters - Stars: individual: R Leo 
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1. Introduction 

R Leo is one of the brightest, apparently largest, and hence best 
studied Mira variable stars. Mira stars are cool, low-mass, pul- 
sating variables that are located on the asymptotic giant branch 
of the Hertzsprung-RusseU diagram and that exhibit a conspic- 
uous mass-loss. Because of to the low temperatures, molecules 
are present in their extended atmospheres, and dust is formed 
at larger distances from the star. High resolution techniques at 
optical/infrared wavelengths, such as speckle interferometry, 
masked aperture techniques, lunar occupation measurements, 
and long-baseline interferometry, allow the measurement of the 
stellar radii, the center-to-limb intensity variations (CLV) in- 
cluding the effects of molecular layers, and the dust envelopes. 
For instance, Mira star diameters and CLVs have been mea- 
sured by, e.g., Di Giacomo (1991"), Haniff et al. (1995), Burns 
et al. 1998 Perrin et al. ( 1999 1, Young et al. (2009), Hofmann et 



al. (l200Tt . Hofmann et al. (2002 1, Woodruff et al ( l2004t . Perrin 
et al. f 2004V Works by, e.g., Danchi et al. fl994'l. Bedding et 
al. (20011, Monnier et al. (2004 1, Ohnaka et al. (2004i also in- 
clude characterizations of circumstellar dust shells. For non- 
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variable giants up to at least spectral type M4, it was shown 
that observed and model-predicted visibility curves do not dif- 
fer significantly from those of uniform disk (UD) profiles up to 
the first minimum, but result in a diameter that is different by 
a few percent at near-infrared wavelengths (see, e.g the discus- 
sions in Wittkowski et al. 2001 , 2004 1. For Mira stars, however, 
it was observed that the visibility function may deviate already 
in its first lobe from a UD curve, indicating a two- or multi- 
component, or Gaussian-shaped intensity profile. Theoretical 
models including the effects of molecular shells have been cre- 
ated that can explain such CLVs. For more detailed discussions 
on these topics, see the reviews by Scholz (1998, 2003), and 
also Scholz i2001> . as well as the references given above. In 
particular, the P and M self-excited dynamic atmosphere model 
series by Hofmann et al. ( 1998 ), Tej et al. ( 2003b i, Ireland et al. 
( 2004a 2004b ), have been designed to match the well-studied 
nearby Miras o Cet and R Leo. These series are based on funda- 
mental mode pulsation models. The comparison of theoretical 
pulsation models with MACHO observations of long-period 
variables in the LMC (Wood et al. 1999 1, pulsation velocities 
derived from Doppler line profiles (Scholz & Wood 2000), as 
well as recent detailed analyses of interferometric radius mea- 
surements (for instance, Mennesson et al. '2002' Woodruff et 
al. 120041 Perrin et al. 2004) strongly indicate that Mira stars 
are fundamental mode pulsators. 
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Table 1. Previous measurements of R Leo's A'-band UD diam- 
eter, together with the visual variability phase at the time of 
observation. 



0^Q[mas] 




Reference 


33.10 ± 1.30 


0.18 


Di Giacomo et al. <199U 


28.18 ±0.05 


0.24 


Perrin et al. (1999) 


30.68 ± 0.05 


1.28 


Perrin et al. 1 19991 


34.00 ± 2.00 


0.44 


Tej et al. 1 1999) 


30.00 ± 0.30 


0.71 


Monnier et al. <2004> 



Earlier diameter estimates of Mira stars which suggested 
overtone pulsation very likely overestimated the true photo- 
spheric diameter due to the effects of molecular layers at larger 
distances from the star (cf. for instance the discussions in Jacob 
& Scholz'2002' Mennesson et al.'2002' van Belle et al.l 2002l 
Wood ruff et al. .2004. Ireland et al. 2004a .2004b. PeiTin et al. 
I2004> . In addition, distances to Mira stars are often not known 
with a high precision since the parallax is typically only a frac- 
tion of the stellar angular diameter. This leads to an additional 
uncertainty when comparing measured Mira star radii to theo- 
retical radius-period relations. Note, however, that the pulsation 
mode determination in Wood et al. ( 1999) depends mainly on 
the period ratios in multiperiodic long-period variables. Since 
period ratios are unaffected by uncertainties in distance or an- 
gular diameter measurements, this mode determination should 
be robust. 

Often, interferometric measurements cannot probe the CLV 
in detail owing to sparse coverage of the wv-plane. In these 
cases, filter-specific uniform disk radii are often transformed 
into physically more meaningful Rosseland or continuum radii 
using available Mira star atmosphere models (e.g. van Belle et 
al. [T9961 12OO2I Hofmann et al. l2002l Boboltz & Wittkowski 
12004 V Detailed tests of Mira star atmosphere models at differ- 
ent stellar variability phases are desirable in order to increase 
the confidence in using them for such purposes. Recently, 
Woodruff et al. ("20041 have shown that their near-infrared K- 
band interferometric observations of o Cet indicate a CLV that 
differs from a UD profile, and that is consistent with CLV pre- 
dictions by the P and M model series mentioned above. 

In this paper, we present a comparison of near-infrared K- 
band VLTI/VINCI interferometric observations of R Leo with 
predictions by the P and M model series. 

R Leo is an oxygen-rich Mira star with spectral type M6e- 
M8IIIe-M9.5e, a period of 3 10 days, a V magnitude of 4.4-1 1 .3 
(Kholopovet al. 1998), and a mass-loss rate of ~ 1 xlO^^Mo/yr 
(Danchi et al. 1994 Knapp et al. 1998 1. Sloan & Price (1998 1 
measured a relatively low dust emission coefficient of 0.23, i.e. 
the ratio of the total emission of the dust to the total emis- 
sion of the star in the mid-infrared. We use a parallax value of 
8.81 + 1.00 mas as given by Whitelock & Feast (EUnOt . which 
is the weighted average of the values by Gatewood (I1992> 
and Pen-yman & ESA (.1997< . Whitelock et al. (l2000t derived 
a mean bolometric magnitude mboi=0.65 with total (peak-to- 
peak) amplitude AOTboi=0.63. Previous measurements of the K- 
band diameter of R Leo are listed in Tabled together with the 
visual phases at the time of the observations. 



Table 2. Details of our observations (date and time of observa- 
tion, spatial frequency, azimuth angle of the projected baseline 
(E of N)), together with the measured squared visibility am- 
plitudes and their errors. The last column denotes the number 
of successfully processed interferograms for each series. The 
effective wavelength for R Leo is ~ 2.22 fim. 
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2. Observations and data reduction 

The R Leo interferometric data were obtained on 1 & 3 April 
2001 (JD=2452003, stellar phase (^vis = 0.08) and on 20 
January 2002 (JD=2452295, 0vis = 1.02) with the ESO 
Very Large Telescope Interferometer (VLTI) equipped with 
the /T-band commissioning instrument VINCI (Kervella et al. 
E003j . The VLTI test siderostats were used on stations EO and 
GO forming an unprojected ground baseline length of 16 m. 
These data were obtained in the early commissioning phase 
of the VLTI, and are publicly available from the ESO archive. 
Figure \l\ shows the visual lightcurve for the cycles covering 
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Fig. 1. Visual lightcurve of RLeo for the cycles covering our 
observations. The data are from the AAVSO ('+' symbols, 
Waagen 2004 1 and AFOEV (diamond symbols, operated by 
CDS) databases. The solid lines denote the dates of maximum 
visual light as given by the AFOEV database. The dashed and 
dashed-dotted lines denote the dates of our observations. Note 
the post-minimum "bump" (see text). 



our observation dates, as obtained from the AAVSO (Waagen 
12004 1 and AFOEV (operated by CDS) databases. Indicated are 
the dates of our observations as well as the dates of maximum 
visual light, as given by the AFOEV database for these cycles. 
It shows a cycle-to-cycle variation from our first to our sec- 
ond epoch of ~ 0.5 mag near the maxima. Note that the post- 
minimum "bump" which is predicted by self-excited pulsation 
models (Hofmann et al.'T998', Ireland et al. I2004al l2004 b i is 
clearly seen in Fig.^ Tab. |2l lists the details of our observa- 
tions together with the obtained calibrated squared visibility 
amplitudes and their errors. Data were obtained as a series of, 
typically, 500 interferograms with scan lengths of 250 //m and 
fringe frequencies of 295 Hz and 415 Hz. Mean coherence fac- 
tors were computed using the VINCI data reduction software, 
version 3.0, as described by Kervella et al. J2004> . employing 
the results based on the wavelets power spectral density. The 
stars Q-Hya (/T-band UD diameter 0^j-, = 9.10 + 0.91 mas, us- 
ing the calibration by Dyck et al. 1996 1, 50ph (0^^ - 9.76 + 
O.lOmas, Borde et al.|,2002i, 9Cen (0^^ = 5.32 + 0.06 mas, 
Borde et al. l2002t, a nd vCet (O^p = 5.18 + 0.51 mas, af- 
ter Dyck et al. [T996t were used as calibration stars. The cal- 
ibration of the visibility values was performed as described in 
Wittkowski et al. (2004 1, using a weighted average of all trans- 
fer function values obtained during the night. R Leo data were 
calibrated using calibration star data obtained with the same 
fringe frequency while the transfer function errors were com- 
puted using all available values. The listed errors of the cali- 
brated squared visibility values include the uncertainties of the 
adopted calibration stars' diameters, the variation of the trans- 
fer function over the night, as well as the scatter of the single 
scans' coherence factors. For the night 3 April 2001, all cal- 
ibration star data were obtained after all RLeo observations. 
This may cause an additional systematic error of the calibrated 



visibility values which is not included in the listed error bars. 
However, it is very unlikely that such an additional calibration 
uncertainty is larger than the error bars shown. All observations 
cover a narrow range of azimuth angles of the projected base- 
line (70-76 deg. E of N), and hence are not sensitive to possible 
asymmetries. The VINCI sensitivity function covers the range 
1.9 to 2.5 fim. The effective wavelength of our RLeo observa- 
tions is Aq ~ 2.22 yum. 

In order to derive effective temperatures from the mea- 
sured angular radius and the bolometric flux (see Sect. |3] be- 
low), we use the mean bolometric magnitude and its amplitude 
given by Whitelock et al. ( i2000il . For the stellar phases 0.08 
and 1.02 of our observations, we derive bolometric fluxes of 
/boi(phase = 0.08) = (1.98 + 0.59) x 10"^ ergcm^^ s"', and 
/boi(phase = 1.02) = (2.05 + 0.61) x 10"^ ergcm"^ s"\ respec- 
tively. Errors of 30% are considered, taking into account cycle- 
to-cycle variations and the uncertainty of the estimate at phases 
0.08 and 1.02. 

3. Results and comparison to models 

Our measured R Leo squared visibility amplitudes are shown in 
Figs.|2land|3] together with a typical model prediction for each 
epoch based on the P and M model series by Hofmann et al. 
(P998I, Tej et al. (2003bt, and Ireland et al. (|2004a 200-®. 
A detailed description of the models and the calculation of 
synthetic visibility values follows below. The visibility curves 
for UD intensity profiles with best fitting UD diameters are 
shown as well. The obtained UD diameters and their formal 
errors for the April 2001 and January 2002 data are ©^^ - 
28.1 + 0.05 mas and 0^^ = 26.2 + 0.01 mas, with reduced xl 
values of 2.47 and 0. 14, respectively. 

The April 2001 data were obtained for a relatively wide 
range of projected baseline lengths (spatial frequencies be- 
tween 31 and 35 cycles/arcsecond) while the azimuth angle of 
the projected baseline changed only marginally, thanks to the 
position of the star on the sky. This allows the study of the CLV 
without the often inherent simultaneous change of baseline an- 
gle and resulting confusion with possible asymmetries. These 
data show that R Leo's CLV at this epoch significantly deviates 
from a UD intensity profile. The January 2002 data cover only 
a much smaller range of spatial frequencies between 34 and 35 
cycles/arcsecond, and do not allow a discrimination of UD and 
more complex intensity profiles. 

The P and M model series are complete self-excited dy- 
namic model atmospheres that were created to match the pa- 
rameters of the best studied Mira variables oCet and RLeo. 
These series are based on fundamental mode pulsation. First- 
overtone pulsation models, such as the O model series of 
Hofmann et al. (1998 1 are not considered since recent results 
indicate that Mira stars are fundamental mode pulsators (see 
Introduction). For the details of the model calculations, in- 
cluding assumptions and approximations, we refer to Hofmann 
et al. (1998 1. The P and M series differ with respect to the 
mass of the so-called "parent star", which is a hypothetical 
non-pulsating giant that has the same mass and luminosity as 
the Mira variable. The geometric pulsation of the Mira occurs 
around the parent's star Rosseland radius R = /^(tross - 1) 



4 



D. Fedele et al.: The A^-band intensity profile of R Leo 



Table 3. Fit results to our data using models from the P and M series with differences of the visual phase between model and 
observations (A^vis = l</'Modei - ^obs.l) up to 0.1 (not taking into account the cycle). Listed are the obtained best-fitting Rosseland 
diameters, and as alternative information the 1.04jum continuum diameters, with their formal errors, the reduced J values, 
the resulting linear Rosseland radii Rrosh and effective temperatures Teft obtained with our adopted values for n and /boi. The 
uncertainties of Rross and Teff are ~ 15% and ~ 10%, respectively. The last columns provide a comparison to the respective model 
values of /Jross and Teg. 
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(e.g. Ireland et al. I2004air2004bt . The parent star of the here 
considered P and M series has solar metallicity, luminosity 
L/Lo=3470, period 332 days, mass M/]VIo=1.0 (P series) and 
1.2 (M), radius «/Rq=241 (P) and 260 (M) (Hofmann et al. 
1199 8 1. The moderately larger radius of the M-series parent star 
leads to a slightly lower effective temperature of the parent star 
and systematically lower phase-dependent effective tempera- 
tures of the pulsating Mira compared to the P series (see the 
discussion in Ireland et al. 2004b i. The M models tend to ex- 
hibit less pronounced cycle-to-cycle variations than the P mod- 
els and to have more compact atmospheres resulting in only 
moderate differences of K-band contamination but larger dif- 
ferences of high-layer absorption features (cf. Hofmann et al. 
[T9981 Ireland et al. 2004a , 2004b.). 

These models are available for 25 phase and cycle combi- 
nations for the P series (see Ireland et al. I2004a> . and for 20 
combinations for the M series (see Ireland et al. l2004b i. The 
models include molecular shells close to continuum-forming 
layers. The model studies by Ireland et al. (2004a, 2004b I show 
that the true continuum CLVs are relatively close to UD pro- 
files and follow a sinusoidal diameter variation as a function 



of photometric stellar phase. For observations with broad near- 
continuum bandpasses in the near-IR that include molecular 
(mainly water) lines, the molecular shells lead to tail-like or 
protrusion-like extensions of the CLV and to deviations from 
the sinusoidal diameter variation (cf. Tej et al. 2003a i. This ef- 
fect due to molecular shells is smallest at pre-maximum phases 
(~ 0.8) of the photometric variability when continuum photons 
are mostly generated in a narrow region just below the strong 
shock emerging from the stellar interior At these phases, effec- 
tive temperatures are highest, and the visibility curves for broad 
bandpasses are expected to be very close to UD curves up to the 
first minimum (see the example in Fig.|^. For maximum and 
post-maximum phases, however, the molecular contamination 
is higher, the CLVs for broad bandpasses are clearly different 
from UD profiles, and these differences can already be detected 
in the first lobe of the visibility function (see the examples in 
Figs.|3and|3below). 

For each model, we used monochromatic tabulated inten- 
sity values at 46 radii between and 5 parent star radii Rp, for 
wavelengths between 1.99 and 2.40 yum in steps of 0.01//m. 
Broad-band synthetic visibility values for the VINCI sensi- 
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tivity function (~ 1.9-2. 5 jum) were computed as described in 
Wittkowski et al. (2004). The VINCI sensitivity curve includes 
the transmission curves of the VINCI /T-filter, the sensitivity 
of the optical fibers, the detector quantum efficiency, and the 
atmospheric transmission. The wings of the VINCI sensitiv- 
ity function between 1.9 and 2.0 jum, as well as between 2.4 
and 2.5 fim are neglected. A dust shell was not included in the 
modeling. The dust emission coefficient of 0.24 (Sloan & Price 
11998 see Sect. is relatively low, and dust shells as optically 
thin as this, are not expected to contribute significantly to K- 
band visibilities (see e.g.. Woodruff et al. 2004J Monnier et al. 
12004 Ohnaka et al. l2004l see also Bedding et al.'200B. Since 
all observations cover only a narrow range of azimuth angles 
of the projected baseline, our observations are not sensitive to 
asymmetries (see above and Tab. As a result, asymmetric 
shapes of the intensity distribution are not considered for the 
modeling. 

We estimate the uncertainty of the assignment of the model 
stellar variability phases relative to the stellar phases at the 
dates of our observations to be about 0. 1 owing to irregulari- 
ties and cycle-to-cycle variations of the light curves of mod- 
eled and real stars (cf. Ireland et al. 2004a, 2004b). As a result, 
fits were made to our observational data with all models that 
have a (model) phase which differs by up to 0. 1 from the stel- 
lar phases at the dates of our observations. For each model, the 
angular diameter is the only free parameter of the fit. In or- 
der to characterize the angular diameter of the fitted CLV, any 
well-defined reference radius of the model CLV can be used, 
such as the Rosseland radius or the 1.04yum continuum radius. 
Physically most meaningful may be a true continuum radius, 
such as the 1.04jum radius, which is not affected by time vari- 
able molecular contamination (see Hofmann et al.Q_998 Jacob 
& Scholz l200l Ireland et al. l2004all2004b l). In the following, 
the Rosseland radius is mainly used as reference quantity, as is 
usual in the literature. It can be transformed into any other ref- 
erence point of the model CLVs with conversion factors derived 
from the model atmospheric structure, as for instance those 
given in Ireland et al. (2004a 2004b). 

Tab. |3l lists for these models the best-fitting Rosseland di- 
ameters, and as alternative information, the 1 .04 yum continuum 
diameters, with their formal errors and the reduced xl values. 
Shown are also the resulting linear Rosseland radii 7?ross and 
the effective temperatures Teff computed with our adopted val- 
ues for the parallax n and the bolometric fluxes /boi (see above). 
The last columns provide a comparison to the respective model 
values of Rross and T^ff. Our derived /?ros,s values have an un- 
certainty of ~ 15%, and the Teff values of ~ 10%. Taking into 
account the uncertainties of our visibility values and of the de- 
rived Rross and Teft values, all considered model predictions are 
consistent with our measured visibility values and with the de- 
rived /?ross and Teff values, to within l-2cr. The correspondence 
of our obtained linear radii with the model radii also indicates 
that the assumption of fundamental mode pulsation is valid. 
For instance, the (piston-driven) first-overtone E model series 
of Bessell et al. ([1996) predict radii (given in Hofmann et al. 
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11998 ) which are larger by a factor of 1 .5 when compared to the 
fundamental mode P series and by a factor of 1 .4 when com- 
pared to the M series. The (self-excited) overtone O series of 



Fig. 2. Measured R Leo squared visibility amplitudes obtained 
in April 2001 (visual stellar phase 0.08), together with the well 
fitting P22 model prediction. For comparison, the UD curve 
is shown as well. The upper panel shows the whole range of 
spatial frequencies while the lower panel shows an enlargement 
of the measured values. 

Hofmann et al. ( 119981 ) have radii larger by factors 2.1 and 1.9, 
respectively. The mean ©ross values, averaged over the fits to 
all considered models, are used as the final angular diameters. 
The average formal error of ©ross is ± 0.03 mas for the April 
2001 data, and + 0.12 mas for the January 2002 data. In addi- 
tion, a calibration error of ±0.1 mas is considered, based on 
calibrations using different calibrators, and different averaging 
of the transfer function over the night. FinaUy, an uncertainty 
arises from the scatter of the fit results using the different mod- 
els, which amounts to +0.39 mas for the April 2001 data, and 
± 0.72 mas for the January 2002 data. Altogether, the final val- 
ues for the Rosseland angular diameters are 

0Ro,ss(April 2001, phase 0.08) = 28.5 + 0.4 mas and 
©Ro,ss(January 2002, phase 1.02) = 26.2 + 0.8 mas. 

The Rosseland angular diameter at the variability phase closer 
to the maximum ( 1 .02) is smaller by ~ 8 % than that at variabil- 
ity phase 0.08. This is consistent with pulsation models (see, 
e.g. Ireland et al. 2004^ I2004b> . Together with the adopted 
values for n and /boi, these angular diameters for April 2001 
and January 2002 correspond to linear radii of 350;^4g R© and 
320+^° Rq, and to effective temperatures of 2930 ± 270 K and 
3080 + 310K, respectively. 

The xi values in Tab. |3] suggest that the models of the 
second cycle of the P series (P2x models) fit our April 2001 
data marginally better than those of the other model cycles. 
However, these differences are not significant, and a better 
baseline coverage would be needed to clearly discriminate be- 
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Fig. 3. Same as Fig.|2 but for the January 2002 (visual stellar 
phase 1 .02) data and the P30 model. 



tween different model cycles and between the P and M model 
series (see the discussion in Ireland et al. I2004al 12004b '). As 
an example, Figs. |3] and |3] show comparisons of the April 

2001 data to the well fitting P22 model, and of the January 

2002 data to the P30 model prediction, respectively. These 
models match the data well at the respective stellar phases, 
in terms of consistency of the CLV and of the Tss and Rrohs 
values. For comparison, best fitting UD curves are shown as 
well. Fig.|3]shows for illustration for our April 2001 data the 
VINCI bandpass-averaged visibility prediction for a typical 
pre-maximum model, P28 (stellar phase 0.83), which is very 
close to a UD profile, and as expected not consistent with our 
post-maximum measurement. 

Fig. shows the corresponding P22, P28, and P30 CLVs. 
Plotted are the monochromatic CLVs as well as the filter- 
averaged CLV. The Rosseland- and the 1.04//m true contin- 
uum radii are indicated, as well as the radius at which the 
filter-averaged CLV decreases by 50%. For the near-UD pre- 
maximum CLV of P28, these different radii coincide because 
there is little molecular contamination. For the maximum and 
post-maximum models P22 and P30 these three different ra- 
dius definitions for each model differ by 15% and 10%, respec- 
tively, owing to a mix of contributions from continuum and 
line-forming layers in our broad bandpass at these phases of 
the pulsation cycle. 

4. Comparison to other observations 

Our R Leo angular diameter values are consistent with earlier 
/T-band measurements listed in Tab^ when taking into account 
possible phase and cycle-to-cycle variations. 
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Fig. 4. For illustration, measured R Leo squared visibility am- 
plitudes obtained in April 2001 (visual stellar phase 0.08), 
as in Fig. |2] but together with a typical pre-maximum (P28) 
model prediction (visual phase 0.83), which is as expected not 
consistent with our post-maximum measurement {xl=2Al). 
The model-predicted CLV of this pre-maximum model is very 
close to a UD model because it includes mainly the CLVs of 
continuum-forming layers without much molecular contamina- 
tion. Consistently, a near-UD /T-band visibility curve was ob- 
served for R Leo by Monnier et al. (.2004 ) at pre-maximum 
phase 0.71. 



Very recently, Perrin et al. (I2004> presented high-precision 
interferometric measurements of R Leo in four relatively nar- 
row sub-filters of the /T-band obtained in November 2001 at 
variability phase 0.80. They compare these measurements to- 
gether with L broad-band observations from 2000 at phase 0.63 
to an empirical model. This model consists of a stellar uniform 
disk and a single thin molecular layer detached from the star 
with different optical depths at each observed wavelength (9 
free parameters), and leads to a good agreement with the data. 
They obtain a stellar (continuum) diameter of 21.88 + 1.7 mas 
at their pre-maximum variability phase. This diameter estimate 
is consistent with our values of the continuum diameter of 
28. 1 + 0.8 mas and 25.9 + 0.9 mas at variability phases 0.08 and 
1.02, considering that the P and M model series predict diame- 
ter differences between phases 0.8 and 1.1 by a factor of about 
1.5 (P08 compared to PI In; P17n to P21n; M08 to Ml In; Ml 8 
to M21n; see Ireland et al. 2004a" 2004bi. It has to be kept in 
mind as well that the modeling approaches are quite different 
and that the assignment of the model phase to the observational 
phase has some uncertainty (see Sect. |3}. In the here adopted 
approach, the fundamental stellar parameters are the only free 
parameters and the CLV in any bandpass is then predicted as 
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Fig. 5. CLV predictions from the (top) P22, (middle) P28, and 
(bottom) P30 model, corresponding to the model visibility 
curves in Figs.|2E and|3 respectively. The P22 and P30 mod- 
els fit well our April 2001 and January 2002 data, respectively. 
The P28 model, a typical pre-maximum model, is as expected 
not consistent with our post-maximum measurements. The thin 
lines denote the monochromatic CLVs while the thick line de- 
notes the CLV averaged over the VINCI sensitivity function. 
Indicated are also (dashed line) the Rosseland radius, (dashed- 
dotted line) the 1.04/im true continuum radius, as well as (dot- 
ted line) the radius at which the filter-averaged CLV drops by 
50%. The mean SiO maser ring radii (for the 43.1GHz and 
42.8 GHz transitions) measured by Cotton et al. ( I2004> close in 
time to our April 2001 data are indicated by the arrows in the 
top panel. 



a function of phase. The measurements of Perrin et al. J2004t 
also illustrate the importance of narrow-band observations for 
more detailed tests of Mira atmosphere models. 

Woodruff et al. (I2004> recently measured the shape of the 
/T-band visibility curve for o Cet at post-maximum stellar phase 
0.13 and found a similar deviation from a UD curve as for 
our post-maximum R Leo observation at phase 0.08, both of 
which are consistent with the predictions of the P model se- 
ries at such stellar phases. This similarity between the RLeo 
and o Cet CLVs at similar phase is remarkable since the shapes 



and amplitudes of their lightcurves are different (Ireland et al. 
I2()()4all2()04bb . 

Monnier et al. f2004^ report that their measured R Leo K- 
band visibility function at pre-maximum phase 0.71 (around 
February 1, 2000) can be well be described by a UD curve. 
This is consistent with the predictions from the P and M model 
series which show near-UD CLVs for broad bandpasses at pre- 
maximum stellar phases (see Sect. |3. Consistently, Perrin et 
al. J1999> measured a R Leo /T-band visibility function which 
is clearly different from a UD curve at post-maximum stellar 
phase ~ 0.28. 

Cotton et al. (2004) have obtained VLBA maps of the 
V = 1,7 = 1-0 (43.1 GHz) and V = 2,7 = 1 - (42.8 
GHz) SiO maser emission toward several Mira stars, including 
RLeo. One of their RLeo epochs, from 29 April 2001, is very 
close in time to our first VLTI epoch (1 & 3 April 2001), with 
a difference of the stellar phase of only 0.09. They derived a 
mean diameter of the SiO maser rings of 51.1 mas and 49.3 mas 
for the V = 1 and v - 2 transitions, respectively (with widths 
of 2.2 mas). Compared to our value for the April 2001 RLeo 
Rosseland angular diameter, the SiO maser spots for the v = 1 
and V - 2 transitions lie at ~ 1.8 and ~ 1.7 Rosseland radii, re- 
spectively. Cotton et al. ( 2004 1 obtained slightly larger mean 
RsiolRi.iijm values of 2.0 and 1.9, respectively, when compar- 
ing to a mean value of the near-infrared /T-band UD diameters 



from Perrin et al. J1999t obtained in April 1996 and March 
1997. Our newly estimated RLeo /?sio/^Ros,s values are also 
in good agreement with recent Rsio/Rross values of ~ 1.7 and 
~ 1.6 for the same v = 1 and v - 2 SiO transitions obtained 
by Boboltz & Wittkowski (2004) for the Mira star S Ori by 
quasi-simultaneous VLTI/VLBA observations. For illustration. 



the mean maser ring radii derived by Cotton et al. ( I2004t are in- 
dicated in the top panel of Fig.|5] together with the well fitting 
P22 near-infrared /T-band model CLV. 
Danchi et al. 



I1994t obtained an inner radius of the dust 
envelope of R Leo of 70 mas by radiative transfer modeling us- 
ing mid-infrared interferometric data taken in 1988 and 1990. 
Compared to our Rosseland angular diameter, this would place 
the inner radius of the dust envelope at roughly 5 Rosseland 
radii, and at roughly three times the distance of the SiO maser 
zone. 



5. Summary and conclusions 

We have obtained near-infrared A'-band interferometric data on 
RLeo for two near-maximum (visual stellar phases 0.08 and 
1.02) epochs. 

For the first epoch, our interferometric data - obtained for 
a sufficiently large range of spatial frequencies - show that 
R Leo's intensity profile deviates significantly from a uniform 
disk profile. We have compared these data to recent self-excited 
dynamic models which were specially created to match the stel- 
lar parameters of o Cet and R Leo. We find that these model 
CLVs for the phases of our observations are consistent with our 
measurements. Furthermore, the model parameters for Rross 
and Teff are consistent with those derived from our observa- 
tions, as well. The models include the effects of close molecular 
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layers, and the correspondence with our measurements show 
that these effects are well described. 

Very recently. Woodruff et al. (2004) found a similar shape 
of the K-hand visibility curve for the Mira star o Cet at post- 
maximum phase 0. 13, and consistency with predictions by the 
P-model series (model P21n). It is remarkable that the two stars 
appear to show similar CLVs which are consistent with the pre- 
dictions by the P model series, while their lightcurves and vari- 
ability amplitudes are different (Ireland et al. ' !2004al I2004bt. 
In accordance with the CLV predictions by the P and M model 
series, Perrin et al. ( 199^ also obtained a R Leo visibility func- 
tion in the /T-band which is very different from a UD curve at 
post-maximum stellar phase 0.28, and Monnier et al. (2004 1 
report a A'-band visibility function which can be well be de- 
scribed by a UD curve at pre-maximum stellar phase 0.71. 
These findings increase our confidence in these dynamic Mira 
star models, which are often used to transform broad-band 
filter-specific UD diameters into more meaningful Rosseland 
or continuum diameters. 

We also obtained high-precision angular Rosseland diam- 
eters for RLeo at the epochs of our observations, and derive 
fundamental parameters, i.e. linear radius and effective tem- 
perature, from these in combination with literature estimates 
of R Leo's distance and bolometric flux. The correspondence 
of our obtained linear radii with model radii of the fundamen- 
tal mode pulsation models used is in agreement with the gen- 
eral recent conclusions that Mira stars pulsate in fundamental 
mode. 

More detailed observations are desirable in the future in or- 
der to better constrain the models. Such observations should 
probe the CLV at a larger range of spatial frequencies. In ad- 
dition, measurements with high spectral resolution in both true 
continuum and certain molecular bands will be useful in or- 
der to separate line-forming and continuum-forming layers. 
Moreover, monitoring of the observed CLVs in time over sev- 
eral cycles with a resolution of ~ 10% of the variability period 
are desirable in order to investigate the strong model-predicted 
CLV variations with variability phase and cycle. 
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